Properties and Mechanisms of Deletions, Insertions, and Substitutions in the Evolutionary History of SARS-CoV-2.

SARS-CoV-2 has accumulated many mutations since its emergence in late 2019. Nucleotide substitutions leading to amino acid replacements constitute the primary material for natural selection. Insertions, deletions, and substitutions appear to be critical for coronavirus's macro- and microevolution. Understanding the molecular mechanisms of mutations in the mutational hotspots (positions, loci with recurrent mutations, and nucleotide context) is important for disentangling roles of mutagenesis and selection. In the SARS-CoV-2 genome, deletions and insertions are frequently associated with repetitive sequences, whereas C>U substitutions are often surrounded by nucleotides resembling the APOBEC mutable motifs. We describe various approaches to mutation spectra analyses, including the context features of RNAs that are likely to be involved in the generation of recurrent mutations. We also discuss the interplay between mutations and natural selection as a complex evolutionary trend. The substantial variability and complexity of pipelines for the reconstruction of mutations and the huge number of genomic sequences are major problems for the analyses of mutations in the SARS-CoV-2 genome. As a solution, we advocate for the development of a centralized database of predicted mutations, which needs to be updated on a regular basis.

Distribution and Functional Analysis of Isocitrate Dehydrogenases across Kinetoplastids.

Isocitrate dehydrogenase is an enzyme converting isocitrate to α-ketoglutarate in the canonical tricarboxylic acid (TCA) cycle. There are three different types of isocitrate dehydrogenase documented in eukaryotes. Our study points out the complex evolutionary history of isocitrate dehydrogenases across kinetoplastids, where the common ancestor of Trypanosomatidae and Bodonidae was equipped with two isoforms of the isocitrate dehydrogenase enzyme: the NADP+-dependent isocitrate dehydrogenase 1 with possibly dual localization in the cytosol and mitochondrion and NADP+-dependent mitochondrial isocitrate dehydrogenase 2. In the extant trypanosomatids, isocitrate dehydrogenase 1 is present only in a few species suggesting that it was lost upon separation of Trypanosoma spp. and replaced by the mainly NADP+-dependent cytosolic isocitrate dehydrogenase 3 of bacterial origin in all the derived lineages. In this study, we experimentally demonstrate that the omnipresent isocitrate dehydrogenase 2 has a dual localization in both mitochondrion and cytosol in at least four species that possess only this isoform. The apparent lack of the NAD+-dependent isocitrate dehydrogenase activity in trypanosomatid mitochondrion provides further support to the existence of the noncanonical TCA cycle across trypanosomatids and the bidirectional activity of isocitrate dehydrogenase 3 when operating with NADP+ cofactor instead of NAD+. This observation can be extended to all 17 species analyzed in this study, except for Leishmania mexicana, which showed only low isocitrate dehydrogenase activity in the cytosol. The variability in isocitrate oxidation capacity among species may reflect the distinct metabolic strategies and needs for reduced cofactors in particular environments.

Blastocrithidia nonstop mitochondrial genome and its expression are remarkably insulated from nuclear codon reassignment.

The canonical stop codons of the nuclear genome of the trypanosomatid Blastocrithidia nonstop are recoded. Here, we investigated the effect of this recoding on the mitochondrial genome and gene expression. Trypanosomatids possess a single mitochondrion and protein-coding transcripts of this genome require RNA editing in order to generate open reading frames of many transcripts encoded as 'cryptogenes'. Small RNAs that can number in the hundreds direct editing and produce a mitochondrial transcriptome of unusual complexity. We find B. nonstop to have a typical trypanosomatid mitochondrial genetic code, which presumably requires the mitochondrion to disable utilization of the two nucleus-encoded suppressor tRNAs, which appear to be imported into the organelle. Alterations of the protein factors responsible for mRNA editing were also documented, but they have likely originated from sources other than B. nonstop nuclear genome recoding. The population of guide RNAs directing editing is minimal, yet virtually all genes for the plethora of known editing factors are still present. Most intriguingly, despite lacking complex I cryptogene guide RNAs, these cryptogene transcripts are stochastically edited to high levels.

Phlebotomine sand fly (Diptera: Phlebotominae) diversity in the foci of cutaneous leishmaniasis in the Surxondaryo Region of Uzbekistan: 50 years on.

In Uzbekistan, the number of reported leishmaniasis cases is rising at the alarming rate. In this work, we studied the phlebotomine sand fly (Diptera: Phlebotominae) diversity in the foci of cutaneous leishmaniasis in the Surxondaryo Region of Uzbekistan and compared it with the data obtained for the same area 50 years ago, when infection prevalence was reportedly low. We found that the implicated vector for zoonotic leishmaniasis, P. papatasi, remained eudominant; the proportion of implicated anthroponotic leishmaniasis vector, P. sergenti, rose significantly from averaged 5.4 to 41.4%; Phlebotomus alexandri, a suspected visceral leishmaniasis vector, was eudominant at two sites, and a second suspected vector for this disease, P. longiductus, was newly recorded in the region. We conclude that the increase in the documented cases of cutaneous leishmaniasis in the Surxondaryo Region of Uzbekistan may be connected to the changes in fauna of sand flies vectoring Leishmania spp.

In silico prediction of the metabolism of Blastocrithidia nonstop, a trypanosomatid with non-canonical genetic code.

BACKGROUND: Almost all extant organisms use the same, so-called canonical, genetic code with departures from it being very rare. Even more exceptional are the instances when a eukaryote with non-canonical code can be easily cultivated and has its whole genome and transcriptome sequenced. This is the case of Blastocrithidia nonstop, a trypanosomatid flagellate that reassigned all three stop codons to encode amino acids. RESULTS: We in silico predicted the metabolism of B. nonstop and compared it with that of the well-studied human parasites Trypanosoma brucei and Leishmania major. The mapped mitochondrial, glycosomal and cytosolic metabolism contains all typical features of these diverse and important parasites. We also provided experimental validation for some of the predicted observations, concerning, specifically presence of glycosomes, cellular respiration, and assembly of the respiratory complexes. CONCLUSIONS: In an unusual comparison of metabolism between a parasitic protist with a massively altered genetic code and its close relatives that rely on a canonical code we showed that the dramatic differences on the level of nucleic acids do not seem to be reflected in the metabolisms. Moreover, although the genome of B. nonstop is extremely AT-rich, we could not find any alterations of its pyrimidine synthesis pathway when compared to other trypanosomatids. Hence, we conclude that the dramatic alteration of the genetic code of B. nonstop has no significant repercussions on the metabolism of this flagellate.

Phylogenetic framework to explore trait evolution in Trypanosomatidae.

The number of sequenced trypanosomatid genomes has reached a critical point so that they are now available for almost all genera and subgenera. Based on this, we inferred a phylogenomic tree and propose it as a framework to study trait evolution together with some examples of how to do it.

Parasites of firebugs in Austria with focus on the "micro"-diversity of the cosmopolitan trypanosomatid Leptomonas pyrrhocoris.

In this work, we investigated parasites of the firebug Pyrrhocoris apterus in Austria and demonstrated that in addition to the extensively studied Leptomonas pyrrhocoris, it can also be infected by Blastocrithidia sp. and by a mermithid, which for the first time has been characterized using molecular methods. This diversity can be explained by the gregarious lifestyle, as well as the coprophagous and cannibalistic behavior of the insect hosts that makes them susceptible to various parasites. In addition, we showed no tight association of the L. pyrrhocoris haplotypes and geographical locations (at least, considering the relatively small scale of locations in Austria) implying that the natural populations of L. pyrrhocoris are mixed due to the mobility of their firebug hosts.

RNA viruses of Crithidia bombi, a parasite of bumblebees.

Leishbuviridae (Bunyavirales) are a diverse monophyletic group of negative-sense single-stranded RNA virus infecting parasitic flagellates of the family Trypanosomatidae. The presence of RNA viruses in trypanosomatids can influence the virulence of the latter. Here, we performed a screening for viruses in Crithidia bombi - a common parasite of important pollinators Bombus spp. (bumblebees) that negatively affects its host in stressful conditions. The majority (8/10) of C. bombi isolates collected in Europe and North America were positive for a virus that we named Crithidia bombi leishbuvirus 1 with high conservation of amino acid sequences between isolates. The results of our comparative phylogenetic analyses of the trypanosomatids and their viruses suggest that the high mobility of bumblebees and frequent coinfections by different strains of C. bombi determine an extensive viral exchange between the latter.

First report of putative Leishmania RNA virus 2 (LRV2) in Leishmania infantum strains from canine and human visceral leishmaniasis cases in the southeast of Brazil.

BACKGROUND: Leishmania RNA virus 1 (LRV1) is commonly found in South American Leishmania parasites belonging to the subgenus Viannia, whereas Leishmania RNA virus 2 (LRV2) was previously thought to be restricted to the Old-World pathogens of the subgenus Leishmania. OBJECTIVES: In this study, we investigated the presence of LRV2 in strains of Leishmania (L.) infantum, the causative agent of visceral leishmaniasis (VL), originating from different hosts, clinical forms, and geographical regions. METHODS: A total of seventy-one isolates were screened for LRV2 using semi-nested reverse transcription-polymerase chain reaction (RT-PCR) targeting the RNA-dependent RNA polymerase (RdRp) gene. FINDINGS: We detected LRV2 in two L. infantum isolates (CUR268 and HP-EMO) from canine and human cases, respectively. MAIN CONCLUSIONS: To the best of our knowledge, this is the first detection of LRV2 in the New World.

Diversity of RNA viruses in the cosmopolitan monoxenous trypanosomatid Leptomonas pyrrhocoris.

BACKGROUND: Trypanosomatids are parasitic flagellates well known because of some representatives infecting humans, domestic animals, and cultural plants. Many trypanosomatid species bear RNA viruses, which, in the case of human pathogens Leishmania spp., influence the course of the disease. One of the close relatives of leishmaniae, Leptomonas pyrrhocoris, has been previously shown to harbor viruses of the groups not documented in other trypanosomatids. At the same time, this species has a worldwide distribution and high prevalence in the natural populations of its cosmopolitan firebug host. It therefore represents an attractive model to study the diversity of RNA viruses. RESULTS: We surveyed 106 axenic cultures of L. pyrrhocoris and found that 64 (60%) of these displayed 2-12 double-stranded RNA fragments. The analysis of next-generation sequencing data revealed four viral groups with seven species, of which up to five were simultaneously detected in a single trypanosomatid isolate. Only two of these species, a tombus-like virus and an Ostravirus, were earlier documented in L. pyrrhocoris. In addition, there were four new species of Leishbuviridae, the family encompassing trypanosomatid-specific viruses, and a new species of Qinviridae, the family previously known only from metatranscriptomes of invertebrates. Currently, this is the only qinvirus with an unambiguously determined host. Our phylogenetic inferences suggest reassortment in the tombus-like virus owing to the interaction of different trypanosomatid strains. Two of the new Leishbuviridae members branch early on the phylogenetic tree of this family and display intermediate stages of genomic segment reduction between insect Phenuiviridae and crown Leishbuviridae. CONCLUSIONS: The unprecedented wide range of viruses in one protist species and the simultaneous presence of up to five viral species in a single Leptomonas pyrrhocoris isolate indicate the uniqueness of this flagellate. This is likely determined by the peculiarity of its firebug host, a highly abundant cosmopolitan species with several habits ensuring wide distribution and profuseness of L. pyrrhocoris, as well as its exposure to a wider spectrum of viruses compared to other trypanosomatids combined with a limited ability to transmit these viruses to its relatives. Thus, L. pyrrhocoris represents a suitable model to study the adoption of new viruses and their relationships with a protist host.

Kinetoplast Genome of Leishmania spp. Is under Strong Purifying Selection.

Instability is an intriguing characteristic of many protist genomes, and trypanosomatids are not an exception in this respect. Some regions of trypanosomatid genomes evolve fast. For instance, the trypanosomatid mitochondrial (kinetoplast) genome consists of fairly conserved maxicircle and minicircle molecules that can, nevertheless, possess high nucleotide substitution rates between closely related strains. Recent experiments have demonstrated that rapid laboratory evolution can result in the non-functionality of multiple genes of kinetoplast genomes due to the accumulation of mutations or loss of critical genomic components. An example of a loss of critical components is the reported loss of entire minicircle classes in Leishmania tarentolae during laboratory cultivation, which results in an inability to generate some correctly encoded genes. In the current work, we estimated the evolutionary rates of mitochondrial and nuclear genome regions of multiple natural Leishmania spp. We analyzed synonymous and non-synonymous substitutions and, rather unexpectedly, found that the coding regions of kinetoplast maxicircles are among the most variable regions of both genomes. In addition, we demonstrate that synonymous substitutions greatly predominate among maxicircle coding regions and that most maxicircle genes show signs of purifying selection. These results imply that maxicircles in natural Leishmania populations remain functional despite their high mutation rate.

Shining the spotlight on the neglected: new high-quality genome assemblies as a gateway to understanding the evolution of Trypanosomatidae.

BACKGROUND: Protists of the family Trypanosomatidae (phylum Euglenozoa) have gained notoriety as parasites affecting humans, domestic animals, and agricultural plants. However, the true extent of the group's diversity spreads far beyond the medically and veterinary relevant species. We address several knowledge gaps in trypanosomatid research by undertaking sequencing, assembly, and analysis of genomes from previously overlooked representatives of this protistan group. RESULTS: We assembled genomes for twenty-one trypanosomatid species, with a primary focus on insect parasites and Trypanosoma spp. parasitizing non-human hosts. The assemblies exhibit sizes consistent with previously sequenced trypanosomatid genomes, ranging from approximately 18 Mb for Obscuromonas modryi to 35 Mb for Crithidia brevicula and Zelonia costaricensis. Despite being the smallest, the genome of O. modryi has the highest content of repetitive elements, contributing nearly half of its total size. Conversely, the highest proportion of unique DNA is found in the genomes of Wallacemonas spp., with repeats accounting for less than 8% of the assembly length. The majority of examined species exhibit varying degrees of aneuploidy, with trisomy being the most frequently observed condition after disomy. CONCLUSIONS: The genome of Obscuromonas modryi represents a very unusual, if not unique, example of evolution driven by two antidromous forces: i) increasing dependence on the host leading to genomic shrinkage and ii) expansion of repeats causing genome enlargement. The observed variation in somy within and between trypanosomatid genera suggests that these flagellates are largely predisposed to aneuploidy and, apparently, exploit it to gain a fitness advantage. High heterogeneity in the genome size, repeat content, and variation in chromosome copy numbers in the newly-sequenced species highlight the remarkable genome plasticity exhibited by trypanosomatid flagellates. These new genome assemblies are a robust foundation for future research on the genetic basis of life cycle changes and adaptation to different hosts in the family Trypanosomatidae.

The 29-nucleotide deletion in SARS-CoV: truncated versions of ORF8 are under purifying selection.

BACKGROUND: Accessory proteins have diverse roles in coronavirus pathobiology. One of them in SARS-CoV (the causative agent of the severe acute respiratory syndrome outbreak in 2002-2003) is encoded by the open reading frame 8 (ORF8). Among the most dramatic genomic changes observed in SARS-CoV isolated from patients during the peak of the pandemic in 2003 was the acquisition of a characteristic 29-nucleotide deletion in ORF8. This deletion cause splitting of ORF8 into two smaller ORFs, namely ORF8a and ORF8b. Functional consequences of this event are not entirely clear. RESULTS: Here, we performed evolutionary analyses of ORF8a and ORF8b genes and documented that in both cases the frequency of synonymous mutations was greater than that of nonsynonymous ones. These results suggest that ORF8a and ORF8b are under purifying selection, thus proteins translated from these ORFs are likely to be functionally important. Comparisons with several other SARS-CoV genes revealed that another accessory gene, ORF7a, has a similar ratio of nonsynonymous to synonymous mutations suggesting that ORF8a, ORF8b, and ORF7a are under similar selection pressure. CONCLUSIONS: Our results for SARS-CoV echo the known excess of deletions in the ORF7a-ORF7b-ORF8 complex of accessory genes in SARS-CoV-2. A high frequency of deletions in this gene complex might reflect recurrent searches in "functional space" of various accessory protein combinations that may eventually produce more advantageous configurations of accessory proteins similar to the fixed deletion in the SARS-CoV ORF8 gene.

Evolution of RNA viruses in trypanosomatids: new insights from the analysis of Sauroleishmania.

RNA viruses play an important role in Leishmania biology and virulence. Their presence was documented in three (out of four) Leishmania subgenera. Sauroleishmania of reptiles remained the only underinvestigated group. In this work, we analyzed the viral occurrence in Sauroleishmania spp. and detected RNA viruses in three out of seven isolates under study. These viruses were of two families-Narnaviridae and Totiviridae. Phylogenetic inferences demonstrated that totiviruses from L. adleri and L. tarentolae group together within a larger cluster of LRV2s, while a narnavirus of L. gymnodactyli appeared as a phylogenetic relative of narnaviruses of Blechomonas spp. Taken together, our work not only expanded the range of trypanosomatids that can host RNA viruses but also shed new light on the evolution and potential routes of viral transmission in these flagellates.

A neo-functionalized homolog of host transmembrane protein controls localization of bacterial endosymbionts in the trypanosomatid Novymonas esmeraldas.

The stability of endosymbiotic associations between eukaryotes and bacteria depends on a reliable mechanism ensuring vertical inheritance of the latter. Here, we demonstrate that a host-encoded protein, located at the interface between the endoplasmic reticulum of the trypanosomatid Novymonas esmeraldas and its endosymbiotic bacterium Ca. Pandoraea novymonadis, regulates such a process. This protein, named TMP18e, is a product of duplication and neo-functionalization of the ubiquitous transmembrane protein 18 (TMEM18). Its expression level is increased at the proliferative stage of the host life cycle correlating with the confinement of bacteria to the nuclear vicinity. This is important for the proper segregation of bacteria into the daughter host cells as evidenced from the TMP18e ablation, which disrupts the nucleus-endosymbiont association and leads to greater variability of bacterial cell numbers, including an elevated proportion of aposymbiotic cells. Thus, we conclude that TMP18e is necessary for the reliable vertical inheritance of endosymbionts.

A Study of Oribatid Mites as Potential Intermediate Hosts of Anoplocephalid Tapeworms of Tatra chamois and Tatra marmots from the Tatra Mountains, Central Europe, and Report of a New Intermediate Host for Andrya cuniculi, the Parasite of Leporidae.

Tatra chamois (Rupicapra rupicapra tatrica (Blahout 1972)) and Tatra marmot (Marmota marmota latirostris (Kratochvíl 1961)) are significant endemic subspecies of the subalpine and alpine ranges of the Tatra Mountains in Central Europe. In four studied localities in the range of their typical biotopes in Slovakia and Poland, we investigated intestinal parasites of Tatra chamois and Tatra marmots, with an emphasis on anoplocephalid tapeworms. We also studied the occurrence, species diversity, and abundance of oribatid mites as intermediate hosts thereof, and the prevalence of cysticercoid larval stages of anoplocephalid tapeworms in collected oribatids using morphological and molecular methods. Coprological analyses revealed the average positivity of Moniezia spp. in chamois faeces at 23.5% and Ctenotaenia marmotae in marmot samples at 71.1%, with significant differences between the localities under study. Morphological analyses determined the presence of cysticercoids in five oribatid species: Ceratozetes gracilis, Edwardzetes edwardsi, Scheloribates laevigatus, Trichoribates novus, and Tectocepheus velatus sarekensis. This is the first record of T. v. sarekensis as an intermediate host of anoplocephalid tapeworms, as well as the first report of Andrya cuniculi occurrence in the territory of the Tatra Mountains, confirmed also by molecular methods.

Genomic analysis of Leishmania turanica strains from different regions of Central Asia.

The evolution in Leishmania is governed by the opposite forces of clonality and sexual reproduction, with vicariance being an important factor. As such, Leishmania spp. populations may be monospecific or mixed. Leishmania turanica in Central Asia is a good model to compare these two types. In most areas, populations of L. turanica are mixed with L. gerbilli and L. major. Notably, co-infection with L. turanica in great gerbils helps L. major to withstand a break in the transmission cycle. Conversely, the populations of L. turanica in Mongolia are monospecific and geographically isolated. In this work, we compare genomes of several well-characterized strains of L. turanica originated from monospecific and mixed populations in Central Asia in order to shed light on genetic factors, which may drive evolution of these parasites in different settings. Our results illustrate that evolutionary differences between mixed and monospecific populations of L. turanica are not dramatic. On the level of large-scale genomic rearrangements, we confirmed that different genomic loci and different types of rearrangements may differentiate strains originated from mixed and monospecific populations, with genome translocations being the most prominent example. Our data suggests that L. turanica has a significantly higher level of chromosomal copy number variation between the strains compared to its sister species L. major with only one supernumerary chromosome. This suggests that L. turanica (in contrast to L. major) is in the active phase of evolutionary adaptation.

Deletions across the SARS-CoV-2 Genome: Molecular Mechanisms and Putative Functional Consequences of Deletions in Accessory Genes.

The analysis of deletions may reveal evolutionary trends and provide new insight into the surprising variability and rapidly spreading capability that SARS-CoV-2 has shown since its emergence. To understand the factors governing genomic stability, it is important to define the molecular mechanisms of deletions in the viral genome. In this work, we performed a statistical analysis of deletions. Specifically, we analyzed correlations between deletions in the SARS-CoV-2 genome and repetitive elements and documented a significant association of deletions with runs of identical (poly-) nucleotides and direct repeats. Our analyses of deletions in the accessory genes of SARS-CoV-2 suggested that there may be a hypervariability in ORF7A and ORF8 that is not associated with repetitive elements. Such recurrent search in a "sequence space" of accessory genes (that might be driven by natural selection) did not yet cause increased viability of the SARS-CoV-2 variants. However, deletions in the accessory genes may ultimately produce new variants that are more successful compared to the viral strains with the conventional architecture of the SARS-CoV-2 accessory genes.